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We present a pedagogic discussion on the role of novel phases of dense 
baryonic matter in “neutron” stars. Qualitative aspects of the physics that 
drives phase transitions and some of its astrophysical consequences are 
discussed. Observable aspects of neutron star structure and early evolution 
of the newly born neutron star are discussed in some detail. 
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1. Introduction 

Despite decades of study, the question - What is the ground state of 
dense baryonic matter? - remains in the realm of theoretical speculation. 
Advances in our understanding of strongly interacting many-particle sys¬ 
tems and QCD suggests numerous candidate states which have competitive 
free energy. Conventional many-body techniques, including those of lattice 
gauge theory and non-relativistic variational and Greens function methods, 
have yet to make definitive statements regarding the phase structure of mat¬ 
ter at densities of relevance to neutron stars. In these lectures we discuss 
a few novel phases of matter and highlight some of their astrophysical con¬ 
sequences. They include (i) baryonic matter with hyperons, (ii) kaon con¬ 
densation and (hi) normal and superconducting quark matter. The choice 
of topics presented is a synthesis of collaborative research the author has 
undertaken in the past four years. In this sense, the scope of the article is 
limited and no attempt is made to present a comprehensive review. For a 
detailed review of assorted topics relating to the interior of neutron stars 
see Ref. [P. 

* Lectures presented at the Zakopane Summer School on Flavor Dynamics, 2002 
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Terrestrial experiments, such as heavy-ion colliders, have been employed 
to probe new phases of matter that are likely to occur at high energy den¬ 
sity. These experiments have proven to be useful in exploring the high 
temperature region of the phase diagram. They have had limited success 
in probing the high-density, low-temperature phases of QCD. Neutron stars 
are best suited to probe this region and thereby provide information that is 
complementary to terrestrial efforts. The first half of the lectures provide a 
lowbrow overview of plausible new phases inside neutron stars. In the latter 
part we discuss how these new phases will affect observable properties of 
neutron stars. In particular, we focus on the following topics to make con¬ 
tact with observables: (i) mass and radius, (ii) the birth and early evolution 
of newly born neutron stars in a supernova explosion, as evidenced by the 
neutrinos they emit. 


2. Nucleonic Matter 

The thermodynamic properties of a many-particle system consisting of 
strongly interacting nucleons is difficult to calculate from hrst principles, 
even if the underlying Hamiltonian were known. Over the past several 
decades there have been numerous attempts to compute the bulk properties 
of nuclear and neutron-rich matter. These include microscopic many-body 
calculations using realistic nucleon-nucleon potentials and phenomenological 
relativistic and non-relativistic mean-field theories. The former approach 
starts with a potential that provides a good description of the measured 
nucleon-nucleon scattering data and uses variational or Greens function 
techniques to obtain the thermodynamic properties of the many-particle 
system. In mean field models the relation to nucleon-nucleon scattering is 
abandoned in favor of a simpler form of the interaction whose parameters 
are determined by fitting the model predictions to empirical properties of 
bulk nuclear matter at nuclear saturation density. 

Walecka’s held theoretical model in which the nucleons interact with 
(j, p and a mesons provides a reasonable description of nuclei and nuclear 
matter |||]. It has been used extensively to study nuclear properties (for a 
review see Ref. ). The model we discuss supplements the original Walecka 
model with self-interactions between mesons and was introduced by Boguta 
and Bodmer Q. The model Lagrangian is given by 

Cn = - QpNl^TN ■ 

+ ^dpad>^a - -U{a)- 

+ \rnlVpy^^ - 


( 2 . 1 ) 
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where 'L'at is the nucleon field operator with tn the nucleon isospin operator. 
Further, = tutv — gaN<y is the nucleon effective mass. This mass is 
reduced from the free nucleon mass miy hy the scalar field cr, taken to have 
rricj = 600 MeV. The scalar self-interaction term is given by 

U{a) = ^ruNigaNcrf + ^{gaNcr)"^ , (2.2) 

where b and c are dimensionless coupling constants. The kinetic energy 
terms for the vector fields corresponding to the uj and p mesons involve 
= dfjVy - duVfj,, and = d^Ry - dyR^ respectively. 

The five coupling constants, g^jN: g^N, gpN, b, and c, are chosen as in 
Ref. to reproduce five empirical properties of nuclear matter at saturation 
density: the saturation density itself is uq = 0.16 fm“^; the binding energy 
per nucleon is 16 MeV; the nuclear compression modulus is 240 MeV; the 
nucleon effective mass at saturation density is 0.78mAr; and the symmetry 
energy^ is 32.5 MeV. 

The model is solved in the mean-field approximation, wherein only the 
time component of the meson fields have nonzero expectation values. The 
symbols cTo,a;o and po below denote expectation values that minimize the 
free energy. In the absence of neutrino trapping, charge-neutral stellar mat¬ 
ter can be characterized by the baryon and electron chemical potentials 
since these are the only conserved charges. Weak interaction equilibrium 
ensures that Pn — Pp = Pe- The free energy density for such a system is 
given by H 

^nuclear ( 

1 f r^Fp \ 

= :^ yj dk (enik) - Pn)+ J dk k^ {ep{k) - pp)j 

+ I - ^1^0 - mlpl) + U{ao) - , (2.3) 

where P{pB,Pe) is the pressure, pB = Pn is the baryon chemical potential 
and 


en{k) = 

sJkP' + - \gpNPo , 

( 2 . 4 ) 

ep{k) = 

JkP' + + gujNi^o + ypNPo , 

( 2 . 5 ) 


are the neutron and proton single particle energies in the mean field approx¬ 
imation. The corresponding Fermi momenta kpn and kpp, which minimize 


^ The energy that depends on the difference in neutron and proton densities. 
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the free energy at fixed baryon and electron chemical potentials, are given 
by solving 


^nikpri) — fJ'n i 
dpikFp) — l-l'p ■ 


( 2 . 6 ) 


Further, the mean field expectation values for the meson fields that minimize 
Eq. (2.3) are given by 


"mlao = g„N {nsikpu) + ns{kFp)) - , 

ml.LVo = QujN {nB{kFn) + nB{kFp)) , ( 2 - 7 ) 

= \gp[nB{kFp) -UBikFn)) , 

where baryon number and scalar densities, ub and Ug respectively, are given 
by 


^ fkpi 


nB{kF) = (^jv^Ar) = — 

■_ ^ Jo 


dk r 


i=n^p 


ns{kF) = 


r^Fi 


r 


m 


i=n^p 


dk k'^ — j= 

0 ^Jk^ + 


N 


( 2 . 8 ) 


m 


He 2 

N 


Note that in Eq. (^) the electron contribution to the free energy has also 
been included. In bulk matter the condition diJlrmcienr/dne = 0, which 
enforces electric-charge neutrality, uniquely determines //e- The magnitude 
and the density dependence of the electron chemical potential is sensitive 
to the value of the nuclear symmetry energy, parametrized in this model as 
the strength of the isovector interaction. 

The model described above provides us with a simple tool to understand 
the various forces at play in the description of charge-neutral nucleonic mat¬ 
ter in weak equilibrium. To study possible phase transitions we will employ 
this model to represent the nuclear phase. Differences between the simple 
model employed here and detailed many-body calculations can, however, 
be important at high densities. This will influence both the location and 
the nature of the phase transitions between the nuclear and more exotic 
phases to be considered below. Nonetheless, the qualitative aspects of the 
phase transitions are generic and are observed even in the most sophisticated 
many-body descriptions of the nuclear phase. 

Fig. H shows the thermodynamic properties of charge-neutral stellar mat¬ 
ter calculated using the mean-field model. The left panel shows the baryon 
chemical potential, electron chemical potential and the energy per particle 
as a function of the baryon density, normalized to the nuclear saturation 
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Fig. 1. The nuclear equation of state. The left panel shows the density dependence 
of the baryon chemical potential, the electron chemical potential and the energy per 
baryon. The right panel shows the relation between energy density and pressure. 


density no = 0.16 fm“^. The right panel shows energy density as a function 
of pressure and is usually referred to as the Equation of State (EoS). It is 
this relation that is required to solve for the structure of the neutron star 
and therefore determines neutron star the mass and radius. 

3. Novel Phases 

With increasing density the chemical potentials for baryon number and 
negative electric-charge increase rapidly due to the repulsive nature of strong 
interactions at short distances. This furnishes energy for the production of 
strange baryons and the condensation of mesons. At even higher densities 
our knowledge of QCD and its asymptotic behavior leads us to expect that 
the quarks inside nucleons will delocalize and form a uniform Fermi sea of 
quarks. These expectations are borne out by model calculations of these 
phases. In what follows we consider a few of these scenarios in detail. 

3.1. Hyperons 

Fig. 2 shows the variation of chemical potential associated with neutral 
and charged baryons. The thick solid curves are predictions of the mean 
field model described in the previous section and the thin lines correspond to 
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the non-interacting Fermi gas model. The upper curve in each of the above- 
mentioned cases corresponds to baryons with negative charge, the middle 
curve to neutral baryons and the lowest curve corresponds to baryons with 
positive charge. The upper and lower horizontal dashed lines indicate the 
vacuum masses of the and the A hyperons respectively. The location 
where the chemical potentials cross the vacuum masses are also indicated 
for the interacting case. This analysis neglects the strong interactions the 
between hyperons and nucleons. If they were on-average attractive, the 
corresponding second-order phase transitions would occur at lower density. 
On the other hand, if they were repulsive, the transitions would be pushed 
to higher densities. Note also that strong interactions between nucleons 
also plays an important role in determining the transition density. For 
the case of non-interacting nucleons the transitions occur at densities that 
are significantly larger. The baryon chemical potential in all models which 
incorporate strong interactions will be larger than the naive prediction of the 
free Fermi gas model at high densities. This generic feature arises because 
strong interactions at short distances are repulsive. 

These are second-order transitions because we are only considering the 
energy cost for introducing infinitesimal hyperon populations. If the hyperon- 
hyperon interactions are strongly attractive, these phase transitions can be¬ 
come first-order. In this case it is energetically favorable to have many 
hyperons appear together rather than for their population to grow infinites¬ 
imally. 

Hyperon-nucleon interactions are poorly constrained due to the lack of 
experimental data. The exception is the A binding energy in hypernuclei Q. 
This information was first employed in mean field models by Glendenning 

f . Hyperons are introduced into the mean-field model discussed in section 
by minimal coupling to the cr, lo and p mesons, and there arise three 
corresponding couplings constants that must be constrained. The binding 
energy of the A provides one constraint which relates the gna and guLj , since 
the A carries no isospin, and additionally the p meson mean field vanishes 
in isospin symmetric matter. In the simplest mean-field model, the hyperon 
couplings to the vector meson couplings are chosen to be similar to 2/3 of 
that of the nucleons. This choice is partially motivated by simple quark 
number and isospin counting rules. The A binding energy is then used to 
determine the strength of the coupling to the a meson. The model further 
assumes that the all hyperons couple to the mesons with the same strength. 
This is likely to be a serious drawback of this simple model. There are 
several direct and indirect indications coming primarily from studies of 
atoms^ that indicate the experiences a repulsive interaction in nuclear 

^ This is an electromagnetic bound state of E“ and a nucleus 
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Fig. 2. Baryon chemical potentials in dense stellar matter computed with (mean- 
field model, thick curves) and without (free Fermi gas model, thin curves) strong 
interactions are shown. In each case, the upper most curve is for negatively-charged 
baryons, the middle curve is for neutral baryon and the lower curve is for positively- 
charged baryons. 


matter [^. 

To summarize, strong interactions between nucleons result in a rapid in¬ 
crease of the baryon chemical potential at high density. The baryon chemical 
potential is typically larger than the lightest hyperon masses at a baryon 
density between 1-3 times nuclear density. Experimental inputs indicate 
that strong interactions between A particles and nucleons are attractive 
and this acts to lower the density for their appearance. The situation with 
the is less clear since experiment suggests that they experience repul¬ 
sive interactions in nuclei, even though the presence of a negative-charge 
chemical potential in neutron stars favors their appearance. 

Hyperons contribute more to the energy density of matter than to the 
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pressure as compared to nucleons. This is easily understood by noting that 
they have larger masses and have smaller Fermi momenta. The change in 
energy density per unit change in pressure is large relative to the nuclear 
case, and this behavior is usually referred to as softening of the EoS. As we 
will discuss in subsequent sections this will act to reduce both the maximum 
mass and the radius of neutron stars. Hyperons can also influence transport 
properties and the thermal evolution of the star since they furnish new 
degrees of freedom that are less degenerate than the nucleons, which can 
more readily participate in dissipative processes. 


3.2. Kaon Condensation 

A large number of electrons are required to ensure charge neutrality in 
dense nuclear matter. The typical electron chemical potential /Ug is about 
100 MeV at nuclear density. With increasing density, fie increases to keep 
pace with the increasing proton number density. The magnitude of this in¬ 
crease depends sensitively on the density dependence of the isovector inter¬ 
action contribution to the nuclear symmetry energy. In Fig. |3| the density 
dependence of the electron chemical potential in the mean-field model is 
shown as the thick black curve. For reference, the electron chemical poten¬ 
tial for the case of non-interacting nucleons is also shown (thin black line). 
Negatively charged bosons whose rest energy in the medium is less than the 
electron chemical potential will condense via a second-order phase transi¬ 
tion. The extent of condensation will be regulated by repulsive interactions 
between bosons in the s-wave channel at low momenta. 

In the hadronic phase the likely candidates for condensation are the tt~ 
and the K~. In vacuum, pions are significantly lighter than kaons, but this 
situation may be reversed in the dense medium due to strong interactions 
between mesons and nucleons. The physical basis for this expectation is that 
the effective theory of meson-nucleon interactions, such as chiral perturba¬ 
tion theory, must incorporate the repulsive s-wave interactions arising due 
to the Pauli principle between the constituent quarks inside nucleons and 
mesons. Mesons containing uoi d quarks will experience repulsion, while on 
the other hand, mesons containing u ov d quarks and s (i.e, kaons) quarks 
will experience attraction in the nuclear medium.^ Experiments with kaonic 
atoms lend strong support to the aforementioned theoretical expectation of 
attraction |^. 

In Fig. ^ the vacuum pion and kaon masses are shown as the dashed 
lines. If the masses do not change in the medium, the figure indicates 
that Tr~ condensation occurs in the vicinity of nuclear density and that K~ 
condensation does not occur for the densities considered. When interactions 

® This simple explanation is due to David Kaplan. 



Reddy-Zakopane02 printed on February 8, 2008 


9 


Meson Condensation 



Fig. 3. Electron chemical potential, the pion and kaon vacuum and in-medium 
effective masses in dense stellar matter 


with medium are included, however, a uniformly charged pion condensate is 
disfavored due to a weak repulsive s-wave interaction. Instead, a spatially 
varying condensate can be favored due to attractive p-wave interactions (for 
a review see Ref. |1R|)- The kaon-nucleon interaction, on the other hand, is 
strongly attractive. In what follows we only discuss kaon condensation for 
this simple reason. 

The idea that kaons could condense in dense nuclear matter was due 
to Kaplan and Nelson |]T|. Using a simplified SUr{3) 0 SUl{3) chiral 
Lagrangian they showed that the K~ could condense at a density about 
three times nuclear density. Subsequently several authors have studied in 
detail the nature and the role of kaon condensation in neutron star matter 
(for a recent review see Ref. p^ ). 

Here, we will employ a simple schematic potential model for kaon- 
nucleon interactions considered in Ref. |13] to illustrate the salient features. 
The scattering length characterizes the low energy kaon-nucleon in¬ 

teraction, and experiment indicates that aj^-^ ~ —0.4 fm. Following the 








10 


Reddy-Zakopane02 printed on February 8, 2008 


authors of Ref. [13| we will assume that the interaction potential can be 
modeled as a square well. The parameters of the potential, i.e, the depth 
Vo = —122 MeV and range R = 0.7 fm, are determined by fitting to low 
energy nucleon-kaon data (for further details see Ref. |^). Given this po¬ 
tential the change in the effective mass of K~ in a low-density medium of 
neutrons can be computed using the Lenz approximation. In this approxi¬ 
mation the attractive potential energy experienced by a kaon at rest can be 
directly related to the scattering length and is given by 


hhenz — 


2-k 

rriKn 


^K~n ) 


(3.1) 


where rriKn is the reduced mass of the neutron-kaon system and is the 
neutron density. The effective mass of the kaon computed using Eq. o is 
shown in the Fig. In this case, the kaon effective mass equals the electron 
chemical potential when ~ 3 — 4 times nuclear density. This corresponds 
to the critical density for kaon condensation. 

At higher densities, the Hartree or mean field approximation is expected 
to be valid. In this case the attractive potential energy of the kaon cannot 
be related directly to on-shell low energy kaon-nucleon scattering data. The 
Hartree potential is given by 


‘-t/i Q 

Ruartree = — Vq R? Un , (3.2) 

where Vq and R are the depth and range of the K~-n potential. Fig. ^ 
indicates that kaons would condense in the Hartree approximation when 

~ 5 — 6 times nuclear density. 

3.3. Normal Quark Matter 

The occurrence of the novel hadronic phases mentioned above depends 
sensitively on the nature of hadronic interactions and their many-body de¬ 
scriptions. In contrast, the asymptotic behavior of QCD, which requires 
that interactions between quarks become weak with increasing momenta, 
provides strong support to the notion that at sufficiently high densities nu¬ 
cleonic degrees of freedom must dissolve to form a nearly free (perturbative) 
gas of quarks. The precise location of this phase transition will depend on 
model descriptions of both the nuclear and quark phases. This is because 
all model studies indicate that the phase transition occurs at rather low 
densities where perturbative methods do not apply. 

The bag model provides a simple description of quark matter and con¬ 
finement. The model was designed to provide a description of the hadron 
mass spectrum. The basic tenants of the model are a non-trivial vacuum 




LI 



Fig. 4. Pressure v/s baryon chemical potential for several quark matter EoS models. 
The pressure in the nuclear mean field model is also shown for reference 


and nearly free quark propagation in spaces (bags) wherein the perturbative 
vacuum has been restored. This restoration costs energy since it requires 
the expulsion of the vacuum condensates. The restoration energy per unit 
volume is called the bag constant and is denoted as B. The model also 
provides a very simple and intuitive description of bulk quark matter. The 
pressure in the bulk quark phase containing up (u), down (d) and strange 
(s) quarks is due to the kinetic energy density of quarks and a negative bag 
pressure. At zero temperature this is given by 


rkpi ry (pjj; j - 

PBagilJ-u: ^^s) — ^i=u,d,s I TTJ ^ 

Jo 


(3.3) 


where 7 = 2spins x 3coior is the degeneracy factor and kpi is the Fermi 
momentum of each quark flavor. The chemical potential for each flavor 

fjii = kp- + mf , where rrii is the corresponding quark mass. In the limit 
of massless quarks and a common chemical potential for all the quarks, the 
pressure in the bag model has the following simple form 




(3.4) 


where /r = /is/3 is the quark chemical potential. The pressure of bulk quark 
matter computed using the Eq. (|3.3|) for B^^‘^ = 150 MeV and B^^'^ = 200 
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MeV are shown in Fig. |^. The phase with the largest pressure is favored 
and the figure shows that the nuclear-quark transition density is very sen¬ 
sitive to the bag constant. For the case when = 150 MeV (the upper 
most curve), three flavor quark matter is the true ground state of matter 


and nuclear matter is a metastable state |lf:]. For larger values of B the 
transition occurs at higher densities. 

Leading order effects of perturbative interactions between quarks can 
also be incorporated in the bag model at high density. This has the effect 
of renormalizing the kinetic term. The pressure in this case is given by 


Pi 


bag,as (/^) — ^^2 


1 + 


2 


TT 


/i -P: 


(3.5) 


where Og = 9^/47r and g is the QCD coupling. At densities of relevance 
to neutron stars, perturbative expansion in ag is not valid. Nonetheless, it 
is still interesting to note that the order o? calculation of the free energy 


predicts a behavior that is similar to that of bag model |^, . Recently, 

Fraga, Pisarski and Schaffner-Bielich recomputed the equation of state of 
massless quark matter to 0{a‘^). They find that the perturbative result is 
well reproduced by the following bag-model-inspired form for the pressure 
over a wide range of densities relevant to neutron stars []T|, 


Pperturb (0 (Q^s ))/^) — . Peff 

47r 


(3.6) 


where = 0.628 and B^^ = 199 MeV for the specific choice for the 
renormalization subtraction point, A = 2/i . The pressure obtained using 
Eq. ( |3.3| ) is also shown in Fig. |^. In the order calculation, the transition 
occurs at even higher density compared to the case when = 200 MeV. 

To summarize, much like in the case of hyperons and kaons, quark matter 
softens the equation of state. The softening is both due to a larger number 
of degrees of freedom in the quark phase and the bag constant which makes 
a negative contribution to the pressure and a positive contribution to the 
energy density. 


3.4- Superconducting Quark Matter 

Since the early work of Bardeen, Cooper and Schrieffer it has been well 
known that degenerate Fermi systems are unstable in the presence of arbi¬ 
trarily weak attractive interactions at the Fermi surface |17|. The instability 
is resolved by the formation of a Bose condensate of Cooper pairs. As is 
well known, for the case of charged fermions, like electrons, this leads to su¬ 
perconductivity. In analogy, the presence of attractive interactions between 
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quarks will lead to pairing and color superconductivity. This was realized 
several decades ago in early work by B. Barrios and S. Frautschi [1^|. Re¬ 
cent realization that the typical superconducting gaps in quark matter are 
larger than those predicted in these early works has generated renewed in¬ 
terest. Model estimates of the gap at densities of relevance to neutron stars 
suggest that A ~ 100 MeV when n = 400 MeV Excellent accounts of 
these recent findings can be in Ref. [^. In what follows, we provide a brief 
introduction to the subject and emphasize aspects that impact neutron star 
phenomenology. 

It is simple to verify that the one-gluon-exchange (OGE) potential be¬ 
tween quarks is attractive in the antisymmetric (color-triplet) channel, such 
as {rib 2 — bir2)jV2. The potential is given by = —2asl3r. Similarly, 
the OGE quark-quark potential in the symmetric (color-sextet) channel is 
repulsive, and the potential is 14^ = a^/Sr. The attraction in the triplet 
channel can result in s-wave pairing between quarks in spin zero and spin 
one channels. Explicit calculations show that the pairing energy, A, is es¬ 
pecially large for the spin-zero case. This type of pairing can only occur 
between unlike flavors of quarks to ensure that the diquark pair (which is a 
boson) has a wavefunction that is symmetric. 

For three massless flavors the condensation pattern that minimizes the 
free energy is called the color-flavor-locked scheme [21|. Pairing occurs be¬ 
tween states that are on opposite sides of the Fermi surface. For massless 
particles the spin-zero state involves either only left and only right-handed 
quarks. The non-zero condensates in this phase are given by 


^ ^abA = A ((5“(5f - 6^6^) (3.7) 

where a,f3 are color indices, a,b are flavor indices, and i,j are spinor. The 
Levi-Civita symbols in spinor space, in flavor space and in color space en¬ 
sures that the pair has spin-zero, and is anti-symmetric in flavor and color 
space respectively. The last equality of this equation explicitly shows that 
the condensation locks color and flavor indices. 

The choice of the colors and flavors that are locked, is of course, arbi¬ 
trary. For one particular choice the pairing scheme is illustrated in Fig. 

In the color-flavor-locked phase (GFL henceforth) all nine quarks participate 
in pairing. Gonsequently, the energy required to create a quark excitation 
involves the breaking of a Cooper pair and the fermion excitation spectrum 
is characterized by an energy gap ~ 2A. Further, gluons acquire a mass via 
the Higgs mechanism by coupling to the colored condensates. The lowest 
energy collective excitations in this phase correspond to Goldstone bosons. 
The number and spectrum of these Goldstone modes can be understood by 

This picture illustration was suggested by George Bertsch. 


4 
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Fig. 5. Pairing scheme in color-flavor-locked quark matter. The different rows 
(columns) correspond to different colors (flavors). 


noting that the condensate breaks baryon number and chiral symmetries. 
Baryon number is broken because the CFL state is not an eigenstate of the 
baryon number operator. This is completely analogous to superconducting 
state of electrons, where the BCS ground state is not an eigenstate of the 
electron number. 


The CFL state breaks chiral symmetry by a novel mechanism Left 
and right-handed condensates are coupled because the color exchange is a 
vector interaction. Flavor rotations on the L condensate must be compen¬ 
sated by a rotation of the R condensate |21], and therefore chiral symmetry 
is spontaneously broken. The octet of flavor Goldstone bosons that result 
will acquire a mass because quark masses explicitly break chiral symmetry. 
The quantum numbers of these pseudo-Goldstone modes maps onto the me¬ 
son octet in vacuum. For this reason they are often commonly referred to 
as the “pions” and “kaons” (or collectively as “mesons”) of the high-density 
phase. Their masses have been computed in earlier work by matching the 
high density chiral effective theory to perturbative QGD. Unlike the vacuum 
case, the square of “meson” masses here is proportional to the product of 
quark masses | |22|| . This results in an inverted hierarchy, where the pions 
are heavier than the kaons. Explicitly, at high density the masses are given 
by 


m^± = a(m„ -k md)ms 
m]^± = a{mu + ms)md 
mj^o = a(md -k ms)mu , 


(3.8) 
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where a = SA^/vr^/i^. In contrast, the baryon number Goldstone boson 
remains massless. This singlet mode is responsible for the superfluid nature 


of this phase [^|. As will be discussed in detail in subsequent sections, this 
massless mode also plays a crucial role in determining the neutrino opacity 
of the CFL phase. 

Pairing at the Fermi surface also contributes to the pressure, which on 
dimensional grounds, we can expect to be ~ For three massless 

flavors, we can supplement the bag model pressure with the contribution 


due to superconductivity. In the CFL phase it is given |2^ 


-PCFL = 


A 


AV^ - B , 


(3.9) 


which is to be compared with Eq. (3.4). 

It is a reasonable approximation to neglect the u and d quark masses 
at densities of relevance to neutron stars where [i ~ 400 MeV. The strange 
quark mass ~ 200 MeV, on the other-hand, cannot be neglected. The 
difference in Fermi momenta between light quarks and the strange quark is 
~ m^/2^. Thus, when A ^ m^/2//, pairing involving strange quarks will be 
suppressed. In the limit of infinite strange quark mass, i.e, in their absence, 
only light quarks pair. This phase is called 2SC (two-flavor superconductor) 
and it is also characterized by pairs that are antisymmetric in flavor ||19|| . 
Antisymmetry in color space excludes one of the three colors from partici¬ 
pating in the condensation. Thus 5C/(3)coior is broken down to 5f7(2)coior- 
Quarks of a particular color and three gluons remain massless. Further, 
unlike the CFL phase, no global symmetries are broken. The absence of the 
massless Coldstone bosons implies the absence of the superfluidity. In this 
sense, the 2SC phase appears quite unremarkable. 

Early attempts to bridge the “gap” between the 2SC phase and the CFL 
phase can be found in Ref. p^ , . They found that the CFL pairing scheme 

is preserved when m<j ^ y/2fiA. For larger mg a first-order transition to the 
2SC phase occurs. The response of the CFL phase to a finite strange quark 
mass in these earlier works ignored the role of the flavor Coldstone modes. 
Bedaque and Schafer showed that the strange quark mass appeared in the 
effective theory for Coldstone bosons in the form of a chemical potential for 
anti-strangeness [26|. Since rriu < my, is lightest Goldstone mode with 
anti-strangeness in the high-density effective theory (see Eq. (^)). Conse¬ 
quently, when the chemical potential, fig ~ m‘^/2^, exceeds the mass of the 

boson they will condense in the ground state. Using the asymptotic for¬ 


mula for the meson masses given in Eq. (3.8), Bedaque and Schafer showed 
that condensation occurs when (m^A^)^/^. The phase with a 
condensate is distinct from the the CFL phase as it breaks additional sym¬ 
metries. However, since the pairing scheme itself remains unaltered we label 
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this phase CFLFC^ In the presence of quark masses and chemical po¬ 
tentials, the CFL phase is symmetric under U{1)y <SiU{1)q rotations, where 

Y is hypercharge or strangeness and Q is the modified electric-charge in the 
CFL phase. In the CFLiF® phase hypercharge or strangeness symmetry is 
spontaneously broken, resulting in the appearance of a massless Goldstone 
boson. 


Much like the strange-quark mass, the presence of an electric-charge 
chemical potential will act as source to condense charged mesons in the CFL 
phase if its strength exceeds the mass of the lightest charged boson. In the 
CFL phase, K~^ and 7r“ are the lowest energy charged-particle excitations. 
(The K~ excitations cost more energy due to the presence of the anti-strange 
chemical potential induced by the strange quark mass.) The phase diagram 
in the presence of an electric charge chemical potential and the strange 
quark mass exhibits a rich structure. For a detailed discussion of these 
novel meson condensed phases and their roles in the birth and evolution of 
neutron stars see Ref. j^]. 

An interesting feature of BCS-like pairing in simple systems with only 
two degrees of freedom is the locking of the Fermi surfaces in momentum 
space. Consider pairing between massless u and d quarks in the presence of 
an electric-charge chemical potential. The charge chemical potential in the 
absence of the pairing interaction will act to change the relative number of u 
and d quarks. In the presence of pairing the state resists this change. With 
increasing electric-charge chemical potential the system exhibits no response 
up to a critical value. When the difference in chemical potentials between 
the particles participating in the pairing exceeds 2A, the system responds 
via a first-order transition to the unpaired normal state. This rigidity was 
initially thought to ensure equal numbers of u, d and s quarks in the CFL 
phase, despite the strange quark mass |^. A detailed analysis, however, 
shows that rigidity in the color-flavor-locked phase takes on a different form 
[29|. It enforces equality between the color and flavor quark number of those 
colors and flavors that are locked by pairing. For the pairing scheme shown 
in Fig. P rigidity manifests in the form of the following equations |^| 


TT'j' — ^11 • 


(3.10) 


The above relations, combined with the additional condition of local color 
neutrality, ensure equal numbers of u, d and s quarks [^]. Therefore, despite 
the strange quark mass, the CFL and the CFLiL^ phases are electrically 
neutral. No electrons are present and there are no massless charged-particle 
excitations. Thus, with regard to electromagnetism, the CFL phase is a 
transparent insulator |Q. 
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4. Nature of the Phase Transition 

The novel phases discussed above could occur either via a first or second- 
order phase transition. In the latter case, the order parameter for the new 
phase is infinitesimal at the critical density and smoothly grows with in¬ 
creasing density. In these situations the electric-charge chemical potential 
is also able to adjust smoothly to maintain charge neutrality in the bulk 
phase. On the other hand, in the case of first-order phase transitions the 
order parameter characterizing the new phase has a finite value at the onset 
of the transition. The electric-charge chemical potential required to ensure 
neutrality changes discontinuously across the transition. In these situations 
the presence of two conserved charges, namely baryon number and elec¬ 
tric charge, allows for the possibility of phase co-existence, i.e., satisfies the 
Gibbs criterion of equal chemical potentials and pressure. 

This phase with two co-existing phases is often referred to as the mixed 
phase. Charge neutrality is enforced globally. The two co-existing phases, 
in our case the nuclear phase and the denser exotic phase, have oppo¬ 
site electric-charges, and carry positive and negative charge respectively. 
The volume fraction occupied by the exotic phase adjusts to ensure overall 
electric-charge neutrality. Denoting the volume fraction of the denser phase 
as X-, global charge neutrality requires 


\ Qexotic (1 X) Qnuclear — 0 • (^-1) 

This was first noted by Glendenning pl| ] in the context of high density phase 
transitions. The phase co-existence occurs over a finite interval in pressure 
due to the presence of the extra degree of freedom, namely electric-charge. 
However, as we discuss next, the extent and existence of such a mixed phase 
depends on surface tension between the two phases. 

To illustrate the physics of first-order phase transitions and the role of 
surface and Goulomb energies in the mixed phase we consider an explicit 
example. The phase transition from nuclear matter to GFL quark matter is 
first-order. The nuclear phase has no strangeness and the bulk quark phase 
has no electrons. The possibility of phase co-existence between these phases 
was investigated in Ref. We highlight some of the main findings here. 
A schematic phase diagram is shown in Fig. |^. The nuclear and GFL phases 
are favored in the lower-left and upper-right shaded regions, respectively. 
With increasing baryon chemical potential, the electron chemical potential 
in the charge-neutral nuclear phase grows, as shown by the solid curve that 
ends at the point labeled A. Nuclear matter is positively charge in the region 
below this curve and negatively charged in the region above. Similarly, the 
GFL phase is electrically neutral when fig = 0, and is negatively (positively) 
charged when fif, > 0 {fj-e < 0). When //g exceeds the mass of the lightest 
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Fig. 6. A schematic phase diagram showing the nuclear and CFL phases in the 
plane of quark chemical potential /r and effective electron chemical potential /ig. 
Isobars are shown as thin solid lines, and neutrality lines for nuclear and CFL 
matter are thick lines. Each phase is negatively charged above its neutrality line, 
and positively charged below it. Continuation onto the unfavored sheet is shown 
by broken lines 


negatively charged Goldstone mode these modes condense (via a second- 
order phase transition) to provide an additional contribution to the electric- 
charge density |Q. A positively charged nuclear phase and a negatively 
charged CFL phase can coexist along the line formed by the intersection 
of the two shaded regions (from A-D). Along this line, which defines the 
mixed phase, both phases have equal pressure and chemical potentials, i.e., 
they satisfy the Gibbs criteria. 

In Fig. 1^ a more detailed version of Fig. the pressure of the bulk nu¬ 
clear, bulk CFL and mixed phases are shown as a function of /x, the quark 
chemical potential. At intermediate values of /x, the mixed phase has larger 
pressure and is therefore favored over both the nuclear and CFL bulk phases. 
The electron chemical potential, /Xe, required to ensure charge neutrality in 
the bulk nuclear phase, grows with ^ as shown. In the mixed phase neu¬ 
trality requires a positively-charged nuclear phase and a negatively-charged 
CFL phase. This is easily accomplished by lowering /Xg from that required 







(MeV) 


Fig. 7. Behavior of the electron chemical potential and the pressure of homogeneous 
neutral nuclear matter, CFL matter and of the mixed phase, all as a function of 
the quark chemical potential /i. Only bulk free energy is included; surface and 
Coulomb energy is neglected. The mixed phase occurs between A and D. The 
vertical line connecting B and C denotes the /x at which the pressures of neutral 
CFL and nuclear matter are equal. This is where a sharp interface may occur. The 
pressure of the mixed phase exceeds that of neutral CFL or neutral nuclear matter 
between A and D |^ . 


to maintain charge neutrality in the nuclear phase (see Fig. ^). The rate at 
which /ie decreases in the mixed phase is shown in the figure and is obtained 
by requiring equal pressures in the both phases at a given /x. 

As noted earlier, in the mixed phase the Coulomb and surface energy 
costs must be met. The results shown in Fig. ^ ignore these corrections. 
In a simple description of the mixed phase one considers a sharp boundary 
between the coexisting phases. A unit cell of the mixed phase is defined as 
the minimum-size region that is electrically neutral. Three different simple 
geometries are considered: spheres, rods and slabs |^. In each of these 
cases, the surface and Coulomb energy cost per unit volume are given by 

^ d X ctqcd 
ro 

= 27r aemfdix) (AQ)^ rg , 


(4.2) 

(4.3) 
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where d is the dimensionality of the structure {d = 1,2, and 3 for slab, 
rod and droplet configurations, respectively), a is the surface tension, and 
AQ = Qnuciear ~ QcFL+kaons is the charge-density contrast between the two 
phases. The other factors appearing in Eqs.( f4.2y4.3| ) are: x, the fraction of 
the denser phase; rg, the radius of the rarer phase (radius of drops or rods 
and half-thickness of slabs); and fd{x), the geometrical factor that arises in 
the calculation of the Coulomb energy which can be written as|32] 


fd{x) = 


d + 2 


2 — d 
d-2 


-I- X 


(4.4) 


The first step in the calculation is to evaluate ro by minimizing the sum of 
and . The result is 


ro = 


d X ctqcd 


47r aemfdix) (A(5)2 


1/3 


(4.5) 


We then use this value of rg in Eqs. 
Coulomb energy cost per unit volume 


to evaluate the surface and 


E^ 


+ E'^ = - (dvr aem d^ fd{x) 


O', 


2/3 

QCD 


(4.6) 


where x equals y when x ^ OA and (1 — x) when 0.5 ^ y ^ 1. The 
dependence x(/x) is computed using Eq. (0)- 

We must now compare this cost to the bulk free energy benefit of the 
mixed phase. The lowest curve in Fig. ^ shows All, the difference between 
the free energy density of the mixed phase (calculated without the surface 
and Coulomb energy cost) and the homogeneous electrically neutral nuclear 
and CEL phases separated by a single sharp interface, whichever of the two is 
lower. The mixed phase has lower bulk free energy, so Afl, plotted in Fig. 
is negative. The remaining curves in Fig. show the sum of the bulk free 
energy difference An and {E^ + the surface and Coulomb energy cost 
of the mixed phase calculated using Eq. ( |4.6| ) for droplets, rods and slabs. 
We employ three different values of (Tqcd since the QCD surface tension 
is, as yet, not very well known. Careful inspection of the figure reveals 
that for any value of ctqcd, the mixed phase is described as a function 
of increasing density by a progression from drops to rods to slabs of CEL 
matter within nuclear matter followed by slabs to rods to drops of nuclear 
matter within CEL matter. This is the same progression of geometries seen 
in the inner crust of a neutron star [ 321 or in the mixed phase at a first-order 
phase transition between nuclear matter and unpaired quark matter [33| or 
hadronic kaon condensate We have also checked that for ctqcd = 10 
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Fig. 8. The free energy difference between the mixed phase and the homogeneous 
neutral nuclear and CFL phases. In the lowest curve, the surface and Coulomb 
energy costs of the mixed phase are neglected, and the mixed phase therefore 
has the lower free energy. Other curves include surface and Coulomb energy for 
different values of ctqcd and different mixed phase geometry. As ctqcd increases, 
the surface and Coulomb price paid by the mixed phase increases. 


and 40 MeV/fm^, with the mixed phase geometry at any x taken to be 
that favored, the sizes of regions of both phases (ro and its suitably defined 
counterpart) are always ^ 6 fm. In general, uniform regions of charge can 
only exist on length scales small compared to the Debye screening length. 
The Debye screening length in the quark and hadronic phase are typically 
of the size 5 to 10 fm |]^ . When the size of the charged regions becomes 
comparable to the Debye screening length, it becomes important to account 
for spatial variations of the charge density. This will influence the surface 
and Coulomb energy estimates presented in Eq. ([4. 61) . A detailed discussion 
of the importance of these finite-size effects is presented in Ref. . 

For any given ctqcD) the mixed phase has lower free energy than ho¬ 
mogeneous neutral CFL or nuclear matter wherever one of the curves in 
Fig. ^ for that ctqcd is negative. We see that some of the mixed phase will 
survive if ctqcd ~ 10 MeV/fm^ while for ctqcd ^ 4,0 MeV/fm^ the mixed 
phase is not favored for any jjL. This means that if the QCD-scale surface 
tension ctqcd ^ 4,0 MeV/fm^, a single sharp interface will be favored. The 
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interface is characterized by a bipolar charge distribution, resulting in an in¬ 
tense electric field which ensures that the electric-charge chemical potential 
is continuous across it (see Ref. [^] for details). 


5. Mass and Radius of Neutron Stars 

The equations that enforce the condition of hydrostatic equilibrium in 
compact stars where general relativistic effects are important are called the 
Tolman-Oppenheimer-Volkov (TOV) equations |36]. They are given by 


dr 

dM{r) 

dr 


-G M{r) e{P) 


^2 g2 


P 

1 + - 
e 


dvrr^P \ 
M{r)c^) 


( 1 - 


2GM{r) 


rc^ 


-1 


= e{P) 


(5.1) 


where P = P{r) and the equation of state specifies e(P), i.e., the energy 
density as a function of the pressure, and M (r) is the total energy enclosed 
within radius r. For a given central pressure, P{r = 0), the above equations 
can be easily integrated out to the surface of the star, where P = 0, to 
obtain the mass and radius of the object. 

A family of stars all described by the same EoS can be obtained by 
parametrically varying the central pressure and repeating the procedure de¬ 
scribed above. Fig. ^ shows the mass-radius curve obtained in this way for 
two different nuclear equation of states. The curve labeled “Mean-Field” 
corresponds to the model described in section ^ and the curve labeled 
APR is a microscopic many-body calculation, by Akmal, Pandharipande 
and Ravenhall, of the EoS using realistic nucleon-nucleon potentials 
The dashed lines labeled R = Rg and R = 1.52Rs (where Rg = 2GM/(? 
is Schwarzchild radius) are model-independent constraints that require all 
stable stars to be to their right. To the left of the line defined by P = the 
objects must be black holes since the event horizon is outside the object, to 
the left of the line defined by P = 1.52Ps the star requires a very stiff EoS 
for which the speed of sound, Cg = dP/de,exceeds the speed of light p8| . 
The last stable orbit around a neutron star is located at a radius P = 3Rg. 
In some accreting neutron star systems called QPO’s (quasi-periodic oscil¬ 
lators) there is observational evidence of the existence of such orbits. This 
implies that they exist outside the physical radius of the star. This con¬ 
strains stars to lie to the left of the line defined by P = 3Rg. Finally, the 
accurate measurement of the neutron star mass in the Hulse-Taylor binary 
system (PSR 1913-1-16) introduces a further constraint which requires that 
the maximum mass of the stars constructed using any model EoS be at least 


as large as 1.44 Mq 
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The Mass-Radius Plot 



Fig. 9. Mass-radius relation for nuclear stars, where Mq is the mass of our sun. 


The radius of the neutron star as inferred from a point far away (com¬ 
pared to Rs) is Roo = R/ \/l — Rs/R where R is the radius of the star, as 
inferred from the TOV equation. An accurate determination of i?oo would 
greatly constrain the allowed region in the mass radius diagram. For exam¬ 
ple, a precise measurement of an i?oo of 10 km constrains all stars to lie on 
the curve labeled Roo = 10 km in Fig. The discovery of a nearby isolated 
neutron star which appears to emit nearly black-body radiation promises 
to provide valuable information about Roo @|- Assuming that the object 
is emitting black-body radiation, a measurement of the photon energy flux 
F and the distance to the object d uniquely determine Roo as 


Roo — d 



(5.2) 


where (Tsb is the Stefan-Boltzmann constant and Too is the temperature of 
the black-body spectrum inferred at infinity. In the case of RXJ 185635- 
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3754, the x-ray spectrum is fit to good approximation by a black-body spec¬ 
trum. The distance to this object d ~ 117 pc However, the optical 

and x-ray data cannot be reconciled with a single temperature black-body 
spectrum Consequently a stringent constraint of the form described 

by the Roo = 10 km curve is not yet possible O, 43|. 


Phase transitions at high density generically result in an increase in 
the energy density of matter at a given pressure. We will loosely refer to 
this as softening. An EoS whose energy density (as a function of pressure) 
is on-average large compared to the nuclear EoS can be considered to be 
soft. A larger energy density increases the gravitation attraction inside the 
compact object. Consequently, it will require higher central pressures to 
achieve hydrostatic equilibrium for a given mass of the compact object. 
However, with increasing pressure, the energy density increases rapidly in a 
soft EoS. The competition between matter pressure and gravity is heavily 
stacked in favor of gravity in a star with a soft EoS. This results in lower 
maximum masses and smaller radii. 


This trend is shown Eig. II- The mass-radius relationship for a hybrid 
compact object containing nuclear matter at low density and pure quark 
matter at high density is shown (black curve labeled NQ). Other exotic 
possibilities like hyperons (NH), kaon condensation (NK) and phase transi¬ 
tions in which a mixed nuclear-quark phase is favored (NQM) have a similar 
effect on the mass-radius curve. Kaon condensation results in stars that are 
relatively compact while the hyperonic and the nuclear-quark mixed phase 
stars favor modestly larger maximum masses and radii. The rough estimate 
of the region of the M-R diagram which can be populated by exotic stars 
given the uncertainties in these models is also indicated by the dashed box. 
Mass and radius are integral quantities and are sensitive to the EoS over a 
wide range of densities. Therefore, in principle one requires a large number 
of different and simultaneous mass-radius measurements to directly infer 
the EoS. However, our knowledge of the nuclear properties at nuclear sat¬ 
uration density and nuclei provide valuable additional information. This is 
encoded in most realistic equations of state and lead to the generic trends 
described above. The accurate measurement of both the mass and radius 
of star will provide a very stringent constraint on the high density EoS. For 
example, the discovery of a compact object with a R ^8 km would provide 
strong support for the existence of novel phases inside neutron stars. On the 
other hand, the accurate measurement of neutron star mass 1.8 Mq 
would support the absence of phase transitions. While neither of these hy¬ 
pothetical measurements will provide conclusive evidence, they will provide 
valuable guidance to the theoretical models of high-density matter. 
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Fig. 10. Mass-radius relation for several exotic stars. The symbols NK, NH and 
NQM refer to locations where we expect the maximum masses of such stars to lie. 
They represent generic trends rather than precise locations. The mean-field EoS is 
employed to represent the nuclear part of the star. 


6. Supernova Neutrinos 


Neutrinos play an important role in stellar evolution. By virtue of their 
weak interactions with matter, neutrinos provide a mechanism for energy 
loss from dense stellar interiors. In neutron stars neutrinos are responsible 
for most of the energy radiated from their birth in a supernova explosion 
through several thousand years of subsequent evolution. In this section we 
present an overview of some of the nuclear/particle physics issues that play 
a role in understanding the rate of propagation and production of neutrinos 
inside neutron stars. The calculation of these rates is of current interest 
since several research groups are embarking on large-scale numerical simu¬ 
lations of supernova and neutron star evolution ^|. Even moderate 
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changes in the nuclear microphysics associated with the weak interaction 
rates at high density can impact macroscopic features that are observable. 
An understanding of the response of the strongly interacting nuclear medium 
to neutrinos and its impact on neutron star evolution promises to provide 
a means of probing the properties of the dense medium itself |47, |^. 

In this section we discuss neutrino interactions in dense matter contain¬ 
ing nucleons and leptons and neutrino interactions in exotic new phases that 
are likely to occur in the dense inner core of the neutron star. We empha¬ 
size the generic aspects of the microphysics that affect the weak interaction 
rates and present supporting qualitative arguments (For a detailed account 
see Ref. 


6.1. Neutrino Interactions in Nucleonic Matter 


It was realized over a decade ago that the effects due to degeneracy 
and strong interactions significantly alter the neutrino mean free paths and 
neutrino emissivities in dense matter |51, 53, 541. However, it is only 


recently that detailed calculations have become available 


57, 58 


61, 62 


64, 


The scattering and absorption reactions that 


contribute to the neutrino opacity are 


+ 73 —> e + B', Vf. + B ^ e'^ + B', 

vx + B ux + B', nx + e~ ^ nx + e~ , (6.1) 


where the scattering reactions are common to all neutrino species and the 
dominant source of opacity for the electron neutrinos is due to the charged- 
current reaction. The important neutrino producing reactions in the neu¬ 
tron star context are 


e + p ^ n + Of., n ^ e + p-\- Vg , 

n +n ^ n + p + e~ + He , n + n ^ n + nvx + Kk (6.2) 


The weak interaction rates for the reactions in Eq. (6.1) and Eq. 
are modified in hot and dense matter because of many in-medium effects. 
The most important of these are: 


(1) Composition: The rate for neutrino processes depends sensitively on the 
composition of the medium, which is itself sensitive to the nature of strong 
interactions. First, the different degeneracies of the different fermions deter¬ 
mine the particle-hole response due to Pauli blocking. For example, a larger 
symmetry energy favors higher proton fractions. This directly impacts the 
weak rates because neutrinos couple differently to different baryonic species 
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and because the Pauli and momentum-conservation restrictions on rates in¬ 
volving neutrons and protons are relaxed. Consequently, the net rates will 
depend on the individual concentrations. 


(2) In-medium dispersion relations'. At high density the single-particle dis¬ 
persion relations are significantly modified from their non-interacting forms 
due to effects of strong interactions. Interacting matter features smaller ef¬ 
fective baryon masses and energy shifts relative to non-interacting matter. 
This in turn affects the weak interaction rates primarily because it modifies 
the density of particle-hole states at the Fermi surface. 


(3) Correlations'. Low-energy neutrinos couple mainly to the long wave¬ 
length density and spin-density fluctuations of the strongly interacting nu¬ 
clear plasma. Repulsive particle-hole interactions and Coulomb interactions 
generally result in a suppression of the weak interaction rates since they in¬ 
crease the energy cost associated with of such fluctuations. On the other 
hand, interactions can also result in low-lying collective excitations to which 
neutrinos can couple. This acts to increase the weak interaction rates at low 
energy. Both of these effects may be calculated using the Random Phase 
Approximation (RPA), in which particle-hole ring diagrams are summed to 

, |64i indicate that 


all orders. Model calculations |51, 53 


, 62 


at high density the neutrino cross sections are suppressed relative to the case 
in which these effects are ignored. In addition, these correlations enhance 
the average energy transfer in neutrino-nucleon collisions. Improvements 
in determining many-body dynamic form factors and assessing the role of 
particle-particle interactions in dense matter at finite temperature are nec¬ 
essary before the full effects of many-body correlations may be ascertained. 


The relative importance of the various effects described above on neu¬ 
trino transport is only beginning to be studied systematically. As a first 
step, we will focus on effects due to modihcations (1) through (3) above. 
To see how this is accomplished we start with a general expression for the 
differential cross section ^ for the processes shown in Eq. o 

1 dV _ Gl Es 

V d^QsdEs ~ ~ 1287r2 ^ 

X (1 - MEs)) Im , (6.3) 

where the incoming neutrino energy is Ei, the outgoing lepton energy is E 3 
and the energy transfer qq = E 3 — Ei. The factor [1 — exp ((—(70 — R 2 + 
F 4 )/kBT)]~^ maintains detailed balance for particle which are in thermal 
equilibrium at temperature T and in chemical equilibrium with chemical 


1 — exp 


-go - [F2 - Fi) 
ksT 
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potentials //2 and /i 4 , respectively. The final state blocking of the outgoing 
lepton is accounted for by the Pauli blocking factor (1 — fsiE^)), where /s 
is the Fermi distribution function. The lepton tensor is given by 

= 8[2k^k^ + {k • q)g^f^ - {k^q^ + q'^k^) ^ ie^^^'^k^q '^], (6.4) 

where is the incoming neutrino four-momentum and is the four- 
momentum transfer. In writing the lepton tensor, we have neglected the 
electron mass terms, since typical electron energies are of the order of a few 
hundred MeV. 

The target-particle retarded polarization tensor is 

= tanh im , (6.5) 


where is the time-ordered (or causal) polarization and is given by 

n„g(g) = -i I ^Tr [T{S 2 {p)r^S,{p + g)r^)]. (6.6) 


The Greens’ functions Si{p) (the index i labels particle species) describe 
the propagation of baryons at finite density and temperature [^]. The 
current operator T^ is 7 ^ for the vector current and 7/^75 for the axial 
current. Effects of strong and electromagnetic correlations may be included 
by calculating the RPA polarization tensor 

j^RPA ^ ^ , ( 6 . 7 ) 


where D denotes the strong interaction matrix (see Ref. |62] for more de¬ 
tails). 

Neutrino mean free path, which is inversely related to the cross section 
per unit volume, calculated in relativistic RPA is shown in Fig. The 
model employed incorporates interactions via a, w, and p exchange as in 
Walecka model 1^, I 


p^ ] . It is supplemented by pion exchange and a repulsive 
contact term, whose strength is parameterized by a constant g' ~ 0.6 to 
account for short-range spin-isospin correlations. The results indicate that 
RPA corrections are most significant in the spin-isospin channel and that 
low-temperature correlation can suppress the cross section by a factor of two 
or more. This is because, in the Fermi-gas limit, the axial-vector response 
of medium makes a contribution that is roughly three times larger than the 
vector response |M]. Quantitative aspects of the suppression depend on 
the details of the model employed. Nonetheless, we note that most model 
studies thus far indicate similar suppression factors. 
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Fig. 11. Neutrino mean free paths in dense stellar matter containing neutron, pro¬ 
tons and electrons calculated in the relativistic RPA {\rpa) and normalized to 
those computed in the mean field or Hartree approximation {\h) ||^ 


6.2. Neutrino Interactions in Novel Phases 

In this section we explore how phase transitions impact the weak in¬ 
teraction rates. As discussed in earlier sections, novel phases of baryonic 
matter are expected to occur at densities accessible in neutron stars. These 
new phases include pion condensation, kaon condensation, hyperons and 
quark matter. An understanding of how these phases might influence neu¬ 
trino propagation and emission is necessary if we are to inquire if these 
phase transitions even occur inside neutron stars. To illustrate how they 
neutrino rates we consider three specific examples of phase transitions: (1) 
generic first-order transitions; (2) superconducting quark matter and (3) 
color-flavor-locked superconducting quark matter. 
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6.2.1. Inhomogeneous Phases: Effects of First-Order Transitions 


First-order phase transitions in neutron stars can result in the forma¬ 
tion of heterogeneous phases in which a positively-charged nuclear phase 
coexists with a negatively-charged new phase which is favored at higher 
densities [^]. This is a generic feature of first-order transitions in matter 
with two conserved charges. In the neutron star context these correspond 
to baryon number and electric-charge. Reddy, Bertsch and Prakash |]66| 
have studied the effects of inhomogeneous phases on z^-matter interactions. 
Based on simple estimates of the surface tension between nuclear matter 
and the exotic phase, typical droplet sizes range from 5 — 15 fm |^|, and 
interdroplet spacings are several times larger. The propagation of neutrinos 
whose wavelength is greater than the typical droplet size and less than the 
interdroplet spacing, i.e., 2 MeV < Ey < 4Q MeV, will be greatly affected 
by the heterogeneity of the mixed phase, as a consequence of the coherent 
scattering of neutrinos from the droplet. 

The Lagrangian that describes the neutral-current coupling of neutrinos 
to the droplet is 

^ -75)z^ (6.8) 


where is the neutral weak current carried by the droplet and Gp = 
1.166x 10“® GeV“^ is the Fermi weak coupling constant. For non-relativistic 
droplets, has only a time-like component, pw{x), where pw is the excess 
weak charge density in the droplet. The total weak charge enclosed in a 
droplet of radius is N\y = d^x pw{x) and the associated form factor 
is 

r'i'd 

F{q) = {\/Nw) / d^x pw{x) sinqx/qx. (6.9) 

Jo 


The differential cross section for neutrinos scattering from an isolated droplet 
is then 

lin 

— + ( 6 . 10 ) 

In the above equation, Ei, is the neutrino energy and 6 is the scattering 
angle. Since the droplets are massive, we consider only elastic scattering for 
which the magnitude of the three-momentum transfer isq = \/2.Ejy(l—cos 9). 

We must embed the droplets in the medium to evaluate the neutrino 
transport parameters. The droplet radius and the inter-droplet spac¬ 
ing are determined by the interplay of surface and Coulomb energies. In 
the Wigner-Seitz approximation, the unit-cell radius is Rw = (3/47rA^£))^/^ 
where the droplet density is Np. Multiple droplet scattering cannot be 
neglected for E^ ^Ij Rw. The effects of other droplets is to cancel scatter¬ 
ing in the forward direction, because the interference is destructive except 
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at exactly zero degrees, where it only produces a change in the index of 
refraction of the medium. 

These effects are usually incorporated by multiplying the differential 
cross section Eq. o by the static form factor of the medium. The 
static form factor, defined in terms of the radial distribution function of 
the droplets, g[r), is 

S{q) = l-^ND J d^rexpiq.f[g{r) — l]. (6-11) 

The droplet correlations, which determine g{r), arise due to the Coulomb 
force and are measured in terms of the dimensionless Coulomb number T = 
Z‘^e^/{SirRy/kT). Due to the long-range character of the Coulomb force, 
the role of screening, and the finite droplet size, g{r) cannot be computed 
analytically. We use a simple form for the radial distribution function g{r) = 
6 {r — Rw)- This choice for g{r) is equivalent to subtracting, from the weak 
charge density pw, a uniform density which has the same total weak charge 
Nyp as the matter in the Wigner-Seitz cell. Thus, effects due to S{q) may 
be incorporated by replacing the form factor F{q) by 

p ) _ 3 sm^R2^-bR;y)co^ 

[qRwr 

The neutrino-droplet differential cross section per unit volume then follows: 


1 da 
V d cos 9 


Nd cos 9)F‘^{q) 

LOTT 


(6.13) 


Note that even for small droplet density Nd, the factor acts to enhance 
the droplet scattering since, as we shall see below Nw ~ 800. 

To quantify the importance of droplets as a source of opacity, we compare 
with the standard scenario in which matter is uniform and composed of 
neutrons. The dominant source of opacity is then due to scattering from 
thermal fluctuations and, in the non-relativistic limit Eq. (|6.3D reduces to 


1 da 
V d cos 9 


3G| El 
lOvr 


[cy (1 -I- cos 0) -I- (3 — cos 0)c^] 


Un 


ksT 

Fpn 


((3.14) 


where cy and ca are respectively the vector and axial coupling constants 
of the neutron, is the neutron number density, Epn = kp^/2Mn is the 
neutron Fermi energy and T is the matter temperature |53]. The transport 
cross sections that are employed in studying the diffusive transport of neu¬ 
trinos in the core of a neutron star are differential cross sections weighted 
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Fig. 12. Neutrino mean free paths as a function of neutrino energy. Solid lines are 
for matter in a mixed phase containing kaons (left panel) and quarks (right panel), 
and dashed curves are for uniform matter. 

by the angular factor (1 — cos 6). The transport mean free path, for 

a given neutrino energy is given by 

Models of first-order phase transitions in dense matter provide the weak 
charge and form factors of the droplets and permit the evaluation of u~ 
droplet scattering contributions to the opacity of the mixed phase. We 
consider two models, namely the first-order kaon condensate and the quark- 
hadron phase transition. The neutrino mean free paths in the mixed phase 
are shown in the left and right panels of Fig. |^, respectively. The results are 
shown for the indicated values of the baryon density ns and temperature 
T, where the model predicts a mixed phase exists. The kaon droplets are 
characterized by radii « 7 fm and inter-droplet spacings Rw ~ 20 fm, and 
enclose a net weak charge N\y rs 700. The quark droplets are characterized 
by r-ci 5 fm, and Rw ~ H fm, and an enclosed excess weak charge 
Nw ~ 850. For comparison, the neutrino mean free path in uniform neutron 
matter at the same ns and T are also shown. The typical neutrino energy 
for neutrino in thermal equilibrium is given by Ei, ~ ttT (the factor vr here 
is just a convenient approximation). For these energies, it is apparent that 
there is a large coherent scattering-induced reduction in the mean free path. 
At much lower energies, the interdroplet correlations tend to screen the weak 
charge of the droplet, and at higher energies the coherence is attenuated by 
the droplet form factor. 


1 da 
V dcosO 


(6.15) 
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The large reduction in neutrino mean free path found here implies that a 
mixed phase will cool significantly slower than homogeneous matter. Conse¬ 
quently, the observable neutrino luminosity at late times might be affected, 
as it is driven by the transport of energy from the deep interior. 


6.2.2. Effects of Quark Superconductivity 

As discussed previously, recent theoretical works |]l8|, W] suggest that 
quarks form Cooper pairs in medium, a natural consequence of attractive in¬ 
teractions destabilizing the Fermi surface. Model calculations, mostly based 
on four-quark effective interactions (Nambu-Jona-Lasinio-like models), pre¬ 
dict the restoration of spontaneously broken chiral symmetry through the 
onset of color superconductivity at relatively low temperatures. They pre¬ 
dict an energy gap of A ~ 100 MeV for a typical quark chemical potential of 
pL ~ 400 MeV. As in BCS theory, the gap will weaken for T > 0, and at some 
critical temperature Tc there is a (second-order) transition to a normal or 
unpaired quark-gluon plasma. During cooling from an initial temperature 
in excess of the formation of a gap in the fermion excitation spectrum 
in quark matter will influence various transport properties of the system. 


Carter and I have studied its influence on the transport of neutrinos |^|. 

The differential neutrino scattering cross section per unit volume in an 
infinite and homogeneous system of relativistic fermions as calculated in 
linear response theory is given by Eq. ( |6.3| ). The medium is characterized 
by the quark polarization tensor In the case of free quarks, each flavor 
contributes a term of the form 


na/3(g) = -iTVc 


d^p 


Tr [5o(p)r„5o(p -b q)Tf 3 ] , 


(6.16) 


where Sq{p) is the free quark propagator. At finite chemical potential and 
temperature the quark propagator is given by 


Soiptaf = 


( A+(p) A {p) \ 
\Po -Ej Pl-Ej) 


{Ppl^ - M7o) • 


(6.17) 


This is written in terms of the particle and anti-particle projection operators 
A'''(p) and A“(p) respectively, where A^(p) = (l±7o7-j5)/2. The excitation 
energies are simply = \p\ — p for massless quarks and = \p\ -\- p for 
massless anti-quarks. The outer trace in Eq. o is over color and sim- 
plihes to 3. The inner trace is over spin, and the Tq are the neutrino-quark 
vertex functions. Specihcally, the vector polarization is computed by choos¬ 
ing (r^jT^) = (70,7/3). The axial and mixed vector-axial polarizations are 
similarly obtained from (ro,r/3) = (7075,7/375) and (ro,r/3) = (70,7/375), 
respectively. 
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The free quark propagators in Eq. ( |6.16 ) are modified in a superconduct¬ 
ing medium. In calculating these effects, we will consider the simplified case 
of QCD with two quark flavors which obey SU(2) 1 0 SU(2) ji flavor symme¬ 
try, given that the light u and d quarks dominate low-energy phenomena. 
Furthermore we will assume that, through some unspecified effective inter¬ 
actions, quarks pair in a manner analogous to the BCS mechanism 0- The 
relevant consequences of this are the restoration of chiral symmetry (hence 
all quarks are approximately massless) and the existence of an energy gap 
at zero temperature, Aq. As discussed earlier in section this super¬ 
conducting phase is called the 2SC phase. The approximate temperature 
dependence of the gap is taken from BCS theory 


A(r) = Aoyi- , (6.18) 

where the critical temperature Tc — 0.57Ao; this relation has been shown 
to hold for perturbative QCD and is thus a reasonable assumption for non- 
perturbative physics. Since the scalar diquark (in the 3 color representation) 
is the most attractive channel, we consider the anomalous propagator 

F{p)abf 9 = {qJaiP)Cl5qgbi-p)) 

A ( A+(p) , A-(p) ^ ^ 

— if^abS^fg^ I 2 £2^ 2 ^ ' (6.19) 

\Po Po ~ ip / 


Here, a,b are color indices, f,g are flavor indices, eabc is the usual anti¬ 
symmetric tensor and we have conventionally chosen 3 to be the condensate 
color. The quasi-particle energy is = yqpf—and for the anti¬ 
particle = -p(IP -|- |u)2 -|- A^. This propagator is also antisymmetric in 
flavor and spin, with C = —^7072 being the charge conjugation operator. 

The color bias of the condensate forces a splitting of the normal quark 
propagator into colors transverse and parallel to the diquark. Quarks of 
color 3, parallel to the condensate in color space, will be unaffected and 
propagate freely. Their response is characterized by the polarization tensor 
given in Eq. ( 6.16| ), but without the trace over color space since it involves 
only quarks of one color. On the other hand, transverse quark colors 1 and 
2 participate in the diquark and thus their quasi-particle propagators are 
given as 


S(f>7/ = ah? + (p.r-mo)- (6.20) 

yt'o sp Po Sp J 

The appearance of an anomalous propagator in the superconducting 
phase indicates that the polarization tensor gets contributions from both 
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the nor mal qu asi-particle propagators Eq. (|6.20 ) and anomalous propaga¬ 
tor Eq. (16.191) . Thus, to order Eq. (|6.16D is replaced with the two 
contributions corresponding to the normal and anomalous diagrams and is 
given by 


n«/3(Q) = [5'o(p)r„S'o(p-F g)r^] 

+ 2Tl [S{p)T^S{p + q)Tp] 

+ 2Tt [F{p)TaF{p +q)Ti3]} . (6.21) 


The S{p)-F{p) mixed terms vanish and the remaining trace is over spin, as 
the color trace has been performed giving the factors of 2. 

For neutrino scattering we must consider vector, axial, and mixed vector- 
axial channels, all summed over flavors. The full polarization, to be used in 
evaluating Eq. (|6.3|), may be written 


Had = E 

f 


( 6 . 22 ) 


The coupling constants for up quarks are Cy ~ k ~ ^ ~ ^ > 

and for down quarks, Cy = ^ sin^ 6w and where sin^ 0\y = 

0.23 is the Weinberg angle. 

The differential cross section, Eq. (|6.3D, and the total cross section are 
obtained by integrating over all neutrino energy transfers and/or angles. 
Results for the neutrino mean free path, A = E/(t, are shown in Fig. |l^ 
as a function of incoming neutrino energy Ei, (for ambient conditions of 
p = 400 MeV and T = 30 MeV). They show the same energy dependence 
found previously for free relativistic and degenerate fermion matter |^l|; 
A oc l/E^ for Ej, ;§> T and A oc 1/E,y for <C T. The results indicate that 
this energy dependence is not modified by the presence of a gap when A ~ T. 
Thus, the primary effect of the superconducting phase is a much larger 
mean free path. This is consistent with the suppression found in the vector- 
longitudinal response function, IIqJj, which dominates the polarization sum 
Eq. (|6.22|) , at ® < Q- 


6.2.3. Neutrino Interactions with Collective Excitations (Goldstone modes) 

The discussion in the preceding section assumed that there were no low 
energy collective excitations to which the neutrinos could couple. As dis¬ 
cussed in Sec. this is true in the two-flavor superconducting (2SC) phase 
of quark matter. For three flavors and when the strange quark mass is neg¬ 
ligible compared to the chemical potential, the ground state is characterized 
by pairing that involves all nine quarks in a pattern that locks flavor and 
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Fig. 13. Left panel: Neutrino mean free path as a function of neutrino energy E^. 
Right panel: Neutrino mean free paths for E^, = ttT as a function of A{T)/T. 
These results are virtually independent of temperature for T ^ 50 MeV. and A 
denote the mean free paths in the superconducting and normal phases respectively. 


color |^|. Naively we can expect significant differences in the weak inter¬ 
action rates between the normal and the CFL phases of quark matter since 
the latter is characterized by a large gap in the quark excitation spectrum. 
However, as detailed in Sec. |3.4| , diquark condensation in the CFL phase 
breaks both baryon number and chiral symmetries. The Goldstone bosons 
that arise as a consequence introduce low-lying collective excitations to the 
otherwise rigid state. Thus, unlike in the normal phase where quark excita¬ 
tions near the Fermi surface provide the dominant contribution to the weak 
interaction rates, in the CFL phase it is the dynamics of the low-energy 
collective states- the Goldstone bosons -that are relevant. Neutrino inter¬ 
actions with Goldstone modes have been recently investigated by Jaikumar 
et al. |TO| and Reddy et ah |n|. 

There are several articles that describe in detail the effective theory for 
Goldstone modes in CFL quark matter [^]. We will not review them here 
except to note that it is possible to parameterize low energy excitations 
about the SU{3) symmetric CFL ground state in terms of the two fields 
B = and S = exp2f(A • tt//^ -|- tj' / /a)- The B fields represent 

Goldstone modes of broken baryon number H. The T, held of broken chi¬ 
ral symmetry, i.e, the pseudo-scalar octet of Goldstone modes vr, and the 
pseudo-Goldstone boson r]'. 

The massless Goldstone boson associated with spontaneous breaking of 
U (1)_B couples to the weak neutral current. This is because the weak isospin 
current contains a havor singlet component. Although neutrinos couple to 
the havor octet of Goldstone bosons, it turns out that the neutrino mean 
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free path is mostly determined by processes involving the massless baryon 
number Goldstone mode |^]. For this reason, we focus our attention on 
these latter processes. The amplitude for processes involving the U{\)b 
Goldstone boson H and the neutrino neutral current is given by |]7T| 




(6.23) 


where = {E, v'^p) is the modified four momentum of the Goldstone boson 
and V = c/\/3 is its velocity |n|. The decay constants for the U{1)b and the 
pseudo-scalar octet of Goldstone modes have also been computed in earlier 
work 1 ^ and are given by fjj = S/U^/Svr^ and = (21 — 81 n 2 )|U^/ 367 r^, 
respectively. 

The neutrino mean free path due to the reaction v ^ Hv can be calcu¬ 
lated analytically and is given by 




(Eu) 


256 

457r 


u(l — u)^(4 -|- v) 
4(1-1- u)2 


Gl 


fn El 


(6.24) 


Neutrinos of all energies can absorb a thermal mode and scatter into 
either a final state neutrino by neutral current processes like u -\- H ^ u and 
V -\- —> 1 / or via the charged-current reaction into a final state electron 
by the process i/g + tt” —> e“. These processes are temperature dependent, 
since they are proportional to the density of Goldstone modes in the initial 
state. The mean free path due to these processes, which we collectively 
refer to as Gerenkov absorption, can be computed. Reactions involving the 
H boson dominate over other Gerenkov absorption processes due to their 
larger population and stronger coupling to the neutral current . For this 
case we hnd the neutrino mean free path is given by 


1 

Kn^niEu) 

F{v,'y) 


F{v,^)GI fjj El 


128 [u (1 -hu)2 
dvr (1 — v) 


52(7) + 


2v 

(1 -u) 


53(7) 


(1 + ^) 

(l-u) 




where 7 = 2vEu/{l — v)kBT and the integrals gn{l) are dehned by 



X 


n 


exp(yx) — 1 


(6.26) 


In contrast to processes involving the emission or absorption of the Gold¬ 
stone modes by neutrinos, the usual scattering process involves the coupling 
of the neutrino current to two mesons. The amplitude for these processes 
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Fig. 14. Neutrino mean free path in a CFL meson plasma as a function of temper¬ 
ature. The neutrino energy E^, = ttT and is characteristic of a thermal neutrino. 
Ass is the total contribution of Cerenkov (absorption and emission) processes to 
the neutrino mean free path and (j) refers collectively to the octet of Goldstone 
modes. 


vanishes for the H meson and is suppressed by the factor p/fn where p is 
the mode momentum for the flavor octet Goldstone modes. Fig. 14 shows 
the contribution of all Goldstone boson-neutrino processes contributing to 
the neutrino mean free path in the CFL phase, including the dominant 
contribution arising from processes involving the massless mode. 

It is interesting to note that the existence of one massless Goldstone 
mode compensates for the large gap in the particle-hole excitation spectrum. 
The contrast between the findings of the previous section where no low 
energy Goldstone modes coupled to neutrinos to those presented here is 
striking. The mean free path in the CFL phase is surprisingly similar to 
that in normal non-superconducting quark phase. 

7. Conclusions 

Neutron stars are excellent laboratories to study possible phase transi¬ 
tions at high baryon density. We have attempted to provide a glimpse of 
the rich phenomena that may arise if such transitions were to occur inside 
these compact objects. The real challenge lies in being able to identify those 
characteristics of the phase transitions that will uniquely affect observable 
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aspects of the neutron star structure and evolution. The structural proper¬ 
ties of the star depend on the EoS of high-density matter. Generically, we 
found that phase transitions result in softening of the high-density EoS. In 
this case, the maximum mass and radius of neutron stars are lowered. 

Neutrino production and propagation in dense matter directly impact 
neutron stars thermal evolution. The presence of novel phases at high den¬ 
sity influences the low-energy weak interaction rates. For specific examples 
considered in this article we found that the rates could be either greatly en¬ 
hanced or reduced. In some instances, the neutrino rates are left unchanged, 
albeit via novel compensating mechanisms. The theory of core collapse su¬ 
pernova evolution, and neutron star cooling, combined with future obser¬ 
vations of supernova neutrinos and multi-wavelength photons from neutron 
stars has the potential to probe the inner core of these remarkable compact 
objects. 
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